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We study specific properties of particles transport by exploring an exact solvable model, a so- 
called comb structure, where diffusive transport of particles leads to subdiffusion. A performance 
of Levy - like process enriches this transport phenomenon. It is shown that an inhomogeneous 
convection flow, as a realization of the Levy-like process, leads to superdiffusion of particles on 
the comb structure. A frontier case of superdiffusion that corresponds to the exponentially fast 
transport is studied and the log-normal distribution is obtained for this case. 
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Fractional kinetics attracts much attention and various 
aspects of thesephenomena are reflected in the recent re- 
views 0, 0, 0> • It should be underlined that while 
subdiffusion described in the framework of the fractional 
equations has numerous applications 0, 0, 0, 0] , there 
are only few examples on fractional superdiffusion equa- 
tions |3 jEy iO | related to Levy like process in dynamical 
chaos 0, 0] and consideration of Levy process for the 
Langevin equations . As it was mentioned in p| "the 
superdiffusion is far from being completely understood" . 
Indeed, it is difficult to imagine the real physical process 
with the infinite variance jump length during the finite 
time (Levy flight). In the Letter we present an example 
of real physical process which is absolutely analogous to 
the Levy flights. In this example we discuss a behavior of 
observable quantities, namely a position of the center of 
mass of traveling contaminant and the form of the packet 
(the tail of a distribution function). The physical mech- 
anism of this Levy-like effect is an inhomogeneous con- 
vection, and the space-time evolution of the contaminant 
in the presence of the inhomogeneous convection flow is 
studied. To emphasize this relation between the inho- 
mogeneous convection and the Levy process, we consider 
the transport phenomenon on a subdiffusive substrate, 
e.g. on a comb structure. 

Our objective is to explain some complications and a 
possible mechanism of superdiffusion on the comb struc- 
ture depicted in Fig. 1. It is an analogue of a subdiffu- 
sive media where subdiffusion has been already observed 
. A comb model is known as a toy model for a porous 
medium used for exploration of low dimensional perco- 
lation clusters ^2 an d electrophoresis process • It 
should be underlined that conditions or changes that are 
necessary to perform in the Liouville equation in order 
to observe superdiffusion in the comb model are impor- 
tant for understanding the nature of this process from 
the general point of view. We show that a performance 
of Levy-like process along, say x-direction, is due to the 
inhomogeneity of the corresponding a;— component of the 
convection current in the Liouville equation. The diffu- 
sion process in such media (modeled by the comb struc- 
ture) are anomalously slow with the subdiffusive mean 
squared displacement of the order of (x 2 (t)) ~ t a , a < 1. 



FIG. 1: The comb structure 



There are external forces leading to convection. In gen- 
eral case, the velocity of the convection flow is space de- 
pendent, i.e. convection is inhomogeneous. The ques- 
tion under investigation is how the observable shape of 
the initial packet changes, when the space-time evolu- 
tion of the packet corresponds to the convection flow. It 
should be underlined that this problem arises in a va- 
riety of applications starting from transport of external 
species (pollution) in water flows through porous geo- 
logical formations C3 , problems of diffusion and re- 
actions in porous catalysts [16j and fractal physiology 
[T7L HH IT^ | . We find the conditions for the inhomoge- 
neous convection velocity, when the observable packet's 
shape corresponds to the normal diffusion with a = 1 for 
the mean squared displacement, or even superdiffusion 
when (x 2 (t)) ~ t a , a > 1. The inhomogeneity will be 
taken in the form of a power law function ~ v s x s . In 
this case, the infinitely long flights take place during a 
finite fixed time. We will show that the Levy process 
takes place when s > 0. Our primary objective is the 
consideration of a frontier case when s = 1. The corre- 
sponding solution of the Liouville equation for the distri- 
bution function of the transport particles is a so-called 
log-normal distribution. 

First, we describe the comb model. A special behavior 
of the diffusion on the comb structure, depicted in Fig. 
1, is that the displacement in the x-direction is possi- 



2 



ble only along the structure axis (x-axis, say at y = 0). 
Therefore the diffusion coefficient as well as a mobility 
is highly inhomogeneous in the y-direction. Namely, the 
diffusion coefficient is D xx — D6(y), while the diffusion 
coefficient in the y-direction along the teeth is constant 
D yy = D. The Liouvillc equation 



dG 
~dt 



div j = 



(1) 



describes a random walk on this structure with the dis- 
tribution function G — G(t,x,y), and the current j = 
(-£><%) |§, -D^). Therefore, it corresponds to the 
following Fokker-Planck equation 



dG 
~dt 



DS(y) 



d 2 G 
dx 2 



D 



d 2 G 



dy 2 







(2) 



with the initial conditions G(0,x, y) = S(x)S(y) and the 
boundary conditions on the infinities G(t, ±oo, ±oo) = 
and the same for the first derivatives respect to x. Per- 
forming the Laplace and the Fourier transforms with re- 
spect to time and the x-coordinate, correspondingly, one 
obtains the solution for the subdiffusion along the struc- 
ture axis |m 



G(t,x,0) 



D 1 ' 2 



dr exp 



Dt 2 



4Dt 



f 



(3) 



The total number of particles on the structure axis de- 
creases with time 



(G) = / G(t,x,0)dx 



1 



(4) 



Therefore this solution of Eq. © describes the subd- 
iffusion when the number of particles is not conserved. 
In what follows consideration, this point will be bear- 
ing in mind poj . Therefore, this solution corresponds 
to the subdiffusion with the second moment or the mean 
squared displacement along the structure axis of the form 



(x 2 (t)) 



(x 2 (t)G(t,x,0)) 
(G(t,x,0)) 



D 



1/2 



(5) 



First, let us consider a transport problem on the comb 
structure due to the constant convection velocity v x = 
VoS(y). Here and in what follows one bears in mind that 
the dimension of v s x s is Dx. The diffusive term in (J2J) is 
omitted, for simplicity. Then, we have 



dG 
~dt 



v 5(y) 



dG 
dx 



D 



d 2 G 



dy 2 



0. 



(6) 



Now applying the Laplace and Fourier transform, this 
equation is solved exactly with the solution in the form 



G(t,x,0) = D l^L eM-Dx 2 /v 2 t) 



(7) 



and G = for x < 0, since the distribution function must 
be positive. Moreover, the is no transport in the negative 
direction due to the initial conditions. It corresponds to 
the normal diffusion with the second moment 



(x 2 (t)) = ^t. 



(8) 



Therefore, to obtain superdiffusion on the comb struc- 
ture, it is not enough to add a fast process e.g. due to 
the constant convection. In what follows we discuss some 
mechanism for the superdiffusion, where we suggest an 
inhomogeneous convection. 

A superdiffusive Levy-like process corresponds to in- 
finite flights during a finite fixed time. Such a process 
could be a realization of an inhomogeneous current in the 
x-direction. For example, the following two-dimensional 
convective current 

j = (v s x s G(t, x, y)S(y), -DdG(t, x, y)/dy) (9) 

after substituting in the Liouville equation gives some 
modification of Eq. © in the form 



dG -v . ^d 2 G „ 
M+V.G6(J,)-D W = 0, 



where s > and 



dG 



V x G = v s x s ^— + sv s x s 1 G. 
dx 



(10) 



(11) 



After the Laplace transform with respect to the time, one 
presents the solution in the form 



G(p,x,y) = /(p,x)exp [-{p/ D) 1 ' 2 ^ 
where we used that 

d 2 



(12) 



exp 



Substituting Eq. O in Eq. {TUJl, wc obtain then the 
following equation for / = f(p,x) 

v 8 x s f + sv s x s - l f + 2[ P DYI 2 f = 5{x) (13) 

with the following boundary conditions on the infinity 
f(p, oo ) = and f{p 7 x) = for x < 0. For s < 1, the 
solution of Eq. I|13|l has the following form 



/ = f exp 



v s x" 



2(Dp) 



V e (l - s) 



(14) 



where Q(x) is a step function: 0(x) = 1 for x > and 
it is zero for x < 0. It corresponds to some kind of 
superdiffusion along the structure axis where all moments 
arc determined by the following distribution function 



G{t,x,0) = 9(x) 



D2X 1 



V 2 (l — s)V 7rf 3 



exp 



Dx 2 



v 2 (l-s) 2 t 



(15) 
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For example, (x 2 (t)} <x t 1 ^ . When s = 0, Eq. (|15|) co- 
incides with Eq. 10. To avoid deficiency at x = where 
the solution is infinite, one needs to introduce a diffusion 
process, by means the second derivative with respect to 
x. In this case, the equations (|10fl . (|ll|) and (|13fl describe 
the space time evolution of the initial profile of particles 
on the asymptotically large scale x S> 1. Therefore, the 
solutions 114|) and (|15|) are the asymptotic ones. When 
s > 1, it is more convenient to look for a solution on 
a class of generalized functions. The case with s > 1 
needs a special care that is deserved a separate consider- 
ation 21]. When s = 1, it is a frontier case, with a strong 
Levy-like process. The so-called log-normal distribution 
could be realized here as well. 

We consider the case of inhomogeneous transport 
along the structure axis due to the inhomogeneous drift 
with s = 1 and inhomogeneous diffusion of the form 
D xx {x,y) = dx 2 S(y). To this end we add the diffusion 
—dx 2 8(y)dG / dx to the x-component of the current in 
Eq. © . The Fokker-Planck equation @ reads now 



where diffusion along the structure axis is 



(16) 



V x G^~dx 2 ^~(2d~v 1 )x^f+vG. (17) 
ox ox 

Again, performing the Laplace transform with respect to 
the time and presenting the solution of Eqs. 116(1 and i|17|) 
in the form of Eq. 1|12|) we obtain the following equation 
for / = f(p,x) 

- dx 2 f" ~ {2d-v 1 )xf + [2[pD] 1 / 2 + Vl )f = S(x) (18) 

with the same boundary conditions as in Eq. J5J. The 
following solution is obtained: 



fiP, x) = 



2[p J D]V2 



(19) 



where r = \[{l-b) + ^f(b+l) 2 + ^(dp) 1 / 2 ], and we de- 
note b = vi/d and d — D jd 2 . The modulus here indicates 
that the equation is invariant with respect to the inver- 
sion x — > — x. We also used here that —xdS(x)/dx = 
5(x) and x 2 d 2 S(x)/dx 2 — 25(x) on a class of constant 
probe functions. This 5-function solution for the inho- 
mogeneous equation can also be obtained by the Fourier 
transform. In this case the specific solution in the Fourier 
space is fk = (2y / pD)~ 1 , whose the inverse Fourier trans- 
form gives the first term in (|19|) . The constant C\ could 
be specified from the initial conditions. It does not carry 
any important information for the large scale asymp- 
totics, and, without loosing of the generality, it could 
be equaled to the unity: C\ = 1. 

To perform the inverse Laplace transform it is conve- 
nient to use the long time asymptotics t — > oo (p — > 0). 
In this case r ~ 1 + a^/p with a = 2d 1 / 2 /(b + 1). It is 



also convenient to present the second term in (|19|) in the 
form of the exponential: \x\~ r = e ^ ln 1^1. Hence the in- 
verse Laplace transform becomes a standard procedure 
|22|. and we obtain the following asymptotic solution: 



G{t,x,0) 




a 2 In' 



2^Wt 



■it 



(20) 



We would like to admit here that |x| 3> 1 and the loga- 
rithmic function is positive. The first term corresponds 
to the smoothly time-decaying pining distribution, while 
the second term, which is the most important for the 
asymptotic transport, corresponds to a some kind of the 
log-normal distribution. It should be underlined that the 
asymptotic solution of Eq. (|20|l is independent of the 
coefficient d, which determines the diffusion along the 
axis structure. It is natural, since an asymptotic space- 
time evolution of the initial distribution is determined by 
the drift component of the current, but not by diffusion. 
It means that one could consider the convective process 
analogous to Eqs. © and 11U|) instead of the more gen- 
eral consideration in the framework of Eqs. (|16fl and 11711 . 
In this case, of course, G(t, x, y) = Q for x < 0. 

Let us calculate the second moment analogously to Eq. 
©. We have then (G) = [1 + i(d + vx)]/VnDt, and the 
second moment is 

(x 2 G) = (8Vi/a 2 )exp[4i/a 2 ]. 

Finely, the displacement is 



(x 2 G) 
(G) 



a 2t/a 2 



(21) 



This exponentially fast spreading is a result of the Levy- 
like process introduced inside the system by means of the 
inhomogeneous displacement current. There is another 
interesting interpretation of the obtained result. It could 
be considered as the Poisson distribution respect to the 
squared logarithm. Namely, if we denote z = In 2 then 
we have from (|20ll 



P{z)dz 



Hence, the average value of z is 



-.e-^z^dz. 



(22) 



zP{z)dz\/[ / P(z)dz] 
n Jo 



it /a 2 . (23) 



It could be understood as the fact that time approaches 
to the infinity with the same rate as the squared loga- 
rithm of the x-coordinate, namely it ~ a 2 In 2 There- 
fore, returning to the log-normal part of the distribution 
function of Eq. H2U|) we have, approximately, that its 
asymptotic behavior is 



G(t»l,|x|>l,y = 0) 



1 



1x1 In" Ixl 



(24) 
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It is simply to see, that the flux is zero on the infinities. 
Finely, one obtains that all even moments of x diverge on 
the large scale asymptotics. The essential difference be- 
tween the distributions for subdiffusion of Eq. © which 
gives for \x\,t > 1, that G(t, \x\,0) oc exp [-(x 4 /t) 1 / 3 ], 
and the superdiffusion of Eq. I|24J) is obvious. Therefore, 
the main mechanism of the transition from subdiffusion 
in Eq. (JSJ) to superdiffusion in Eq. (|21|l on the comb 
structure is the inhomogeneous convection. 

In conclusion, we would like to classify a possible re- 
sponse of a substrate with anomalous transport. As it 
is shown here it depends on the external forcing accord- 
ing to the power rate of the inhomogeneous convection 
j x (x,t) — v s x s 8{y)G{t,x,Q). As it follows from the so- 
lution (|15|l , when s < it is subdiffusion 0, 153] • As we 
shown here, when s > it is superdiffusion, Eq, i|15[l . The 
homogeneous convection with s — corresponds formally 
to the normal diffusion as in Eq. JSJl, but the effective 
diffusion coefficient Vq/D is determined by the external 



forcing vq. The frontier case with ,s = 1 has two features. 
The first one corresponds to the log-normal distribution 
of transport particles, where one deals not with a sum of 
independent random variables but with their multiplica- 
tion [8j- As the result the exponentially fast transport 
takes place (see Eq. (|21|) ). The second feature is that 
the solution satisfies to the natural boundary conditions, 
when the flux of transporting particles equals zero on 
the infinities for any time. When s > 1, it corresponds 
to superdiffusion [21j, as well. 

The final remark is that the relaxation along the 
structure axis is "inevitable" process due to the item 
DS(y)dG/dx. Therefore, anomalous transport is due 
to the convection flow only on the asymptotically large 
scale x 3> 1 where the diffusion process could be omitted. 
This point is studied in detain in |2l| . where we shown 
that this approximation corresponds to Liouville-Green 
asymptotic solution |24| . 
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